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Topological semimetals are gapless states of matter which have robust surface states and inter¬ 
esting electromagnetic responses. In this paper, we consider the electromagnetic response of gapless 
phases in 3 -I- 1-dimensions with line nodes. We show through a layering approach that an intrin¬ 
sic 2-form Bfiv emerges in the effective response field theory that is determined by the geometry 
and energy-embedding of the nodal lines. This 2-form is shown to be simply related to the charge 
polarization and orbital magnetization of the sample. We conclude by discussing the relevance for 
recently proposed materials and heterostructures with line-node fermi-surfaces. 


Topological insulators (TIs) have been of great inter¬ 
est in recent years after their theoretical proposal and 
experimental discovery in the past decade. Their elec¬ 
tronic properties led to a wide search for novel topolog¬ 
ical band structures in many materials[l, 2]. TIs are 
characterized by a gapped bulk and protected bound¬ 
ary modes that are robust in the presence of disorder. 
They also exhibit quantized properties in their electro¬ 
magnetic (EM) response [3, 4]. A classification of non¬ 
interacting fermionic states protected by discrete time- 
reversal (T), charge-conjugation (C), and chiral symme¬ 
tries has been worked out in Refs. 4-6. This has fur¬ 
ther been expanded on in recent years to include trans¬ 
lation, reflection, and rotation symmetries of crystalline 
systems [7-22]. These theoretical advances have been ac¬ 
companied by experimental discoveries of several TIs in 
various symmetry classes. The 3D T-invariant strong 
TI (e.g., BiSb[23], 612803 [24-26]), the 2D quantum spin 
Hall insulator (e.g., CdTe/HgTe quantum wells [27-29]), 
the 2D quantum anomalous Hall (Chern) insulator (e.g., 
Cr-doped (Bi,Sb) 2 Te 3 [ 3 , 30]), and a 3D topological crys¬ 
talline insulator (PbSnTe)[31, 32]. 

A defining characteristic of TIs is a gapped bulk, but 
one can also ask if there are gapless states of matter 
which harbor protected boundary modes and have un¬ 
usual EM responses and transport properties. This ques¬ 
tion has been asked, and answered in the affirmative 
with the discovery of topological semimetals (TSMs). 
The most studied TSMs all have point-like Fermi sur¬ 
faces, e.g., 2D Dirac semi-metals/graphene [33], 3D Weyl 
semimetals[34, 35], and 3D Dirac semimetals [36-41]. 
In recent work, we proposed a unifying structure to 
understand TSMs with point-like Fermi-surfaces, from 
which one can straightforwardly determine the quasi- 
topological EM responses [42], and which expands on pre¬ 
vious work[43-49]. The main perspective which helps us 
understand these TSMs are models produced by a layer¬ 
ing construction. Generically, a TI phase in d spatial di¬ 
mensions can be layered/stacked into d -|-1 dimensions by 
introducing “trivial” tunneling between the layers, i.e., 
tunneling that does not immediately generate a d -|- 1 - 
dimensional strong topological phase. As the tunneling 


coefficient is increased, we showed that generically the 
material will transition from a weak topological insulator 
phase, which is formed in the decoupled limit, to a trivial 
insulating state with an intervening semimetallic gapless 
phase with point-nodes. To ensure the stability of the 
gapless phase additional symmetries are often required. 
For example, in the case of the 2D Dirac semimetal, the 
stability of the intermediate TSM phase relies on the 
presence of a composite spatial and anti-unitary sym¬ 
metry, i.e., TT where I is inversion symmetry. 

The EM response of point-node TSMs is generally 
characterized by an intrinsic 1 -form b = (bo, hi) which is 
related to the locations of the nodes in momentum and 
energy space[42, 44-46]. This quantity is analogous to 
the weak invariant P = (where Vi are integers, 

and Gi are reciprocal lattice vectors) that exists in the 
related, fully-gapped weak TI phase, but it additionally 
contains an extra time component bo due to the energy 
difference between the nodes. The time-component is 
ill-defined in the gapped phase unless, e.g., the system 
is subjected to a periodic driving frequency. The actual 
dependence of the EM response on b depends on the type 
of point-node semimetal, and can generate a wide variety 
of effects in 2D and 3D TSMs. 

While a TSM with point-nodes arises from coupling d- 
dimensional topological phases into a, d + 1 -dimensional 
system, we can extend this idea by layering a d- 
dimensional topological phase into d + 2 dimensions. 
When the d-dimensional elements are decoupled, the 
d + 2 -dimensional system will be in a secondary weak 
topological phase characterized by an anti-symmetric 
tensor/2-form invariant r'ij[50, 51]. When the lower di¬ 
mensional topological phases are coupled with strong- 
enough “trivial” hopping then they will produce line- 
node Fermi surfaces (FLs). We mentioned in Ref. 42 
that in these gapless phases we expect the EM response 
to be characterized by an analogous 2-form 5^,^. Hence, 
the goal of this article is to conclusively show that the 
effective response action is given by 

S[A,B] = Y^J (1) 
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where is related to the magnetization and polar¬ 
ization via eSoi = and eBij = Air'^eijkP^ for 

i = x,y, z. The 2-form is an intrinsic property of 

line-node semimetals determined by the geometry of the 
nodal submanifolds, and is the analogue of a secondary 
weak invariant, though for a gapless phase. It also in¬ 
cludes components where fi, v are in the time direction, 
which are not available for a time-independent gapped 
system. 

To summarize, in this article, we consider the EM re¬ 
sponse of 3D TSMs with non-degenerate line-like Fermi 
surfaces which we will term “LTSMs.” We note that a 
gapped version of this stacking construction has also been 
considered in [52], while a superconducting version of 
this, including line nodes, has been considered in [53, 54]. 
Additionally, in very recent work, several proposals for 
materials that realize line-node TSM states have ap¬ 
peared which utilize magnetic heterostructures [55, 56], 
carbon allotropes [57, 58], and inversion symmetric 
Cu 3 PdN[ 59 , 60]. After our primary discussion, we com¬ 
ment on how our analysis could be used to generate an 
EM response in these systems, including systems with 
nominally spin-degenerate FLs. 

To aid our discussion it will be helpful to consider an 
explicit model. Let us take the 3D Bloch Hamiltonian 

H^{k) = suvkxCj'^ + — m — coskx — ( 2 ) 

(3 cosky — 'y cos kz)<J^, 

which has inversion X — and time reversal T = 
symmetries, where cr“ represent two (non-spin) degrees of 
freedom, and the lattice constant a = 1. When (3 = j = 
0, this model reduces to decoupled ID insulators aligned 
parallel to the cc-direction. Since each ID wire is inver¬ 
sion symmetric, their polarization will be quantized (and 
all equal). In their topological phase, the polarization of 
a single wire will be Px{ky,kz) = ^ ^ PT[Ax(k)\dkx = 
e/2mode[ll, 61-64], where Ai{kx) is the adiabatic con¬ 
nection matrix Af‘{k) = —i{ua,k\-;^\ub,k), where a, b run 
over the occupied bands. If each insulator was instead 
in a trivial state, we would have Px(ky,kz) = Omode. 
The total polarization is just the sum over all the decou¬ 
pled wires. In addition to the bulk topological proper¬ 
ties, the ID TIs have degenerate mid-gap modes local¬ 
ized at opposite ends of the system, the filling of which 
determines the bound surface charge. To unambigu¬ 
ously determine the sign of the bulk polarization, and 
hence the sign of the surface charge, one must break the 
degeneracy by adding an infinitesimal (inversion) sym¬ 
metry breaking mass, e.g., and take the limit as 

TOx —t 0. Hence, j3 = ^ = Q implies a secondary weak 
TI state protected by inversion symmetry and the EM 
response is given by Eq. 1, but for the special case when 

B = sgcimi{\Gy AGz) => Px = sgnmi^^. 

Now let us tune away from the decoupled limit. When 
7 = 0 and (3 is increased such that the bulk gap closes, we 


get layers of 2D Dirac semimetals with two gapless points 
in the Brillouin zone (BZ) at = ± cos“^ each 

value of kz- This system is in a non-generic LTSM phase 
with two straight FLs that traverse the entire BZ. The 
nodes are locally stable in the BZ as long as the compos¬ 
ite TT-symmetry is preserved. In fact, in terms of the 
flux of the adiabatic connection, we note that the Dirac 
nodes carry a Berry flux of ± 7 r, and hence are stabilized 
from the formation of an energy gap since doing so would 
smoothly spread the flux around the gapped degeneracy 
point, which is forbidden by TI symmetry [2, 42]. The 
FLs carry a helicity y, which along with sgn tox, indicates 
the sign of the Berry flux. The response theory of this 
simple case can be directly determined from the results 
of Ref. 42. Indeed, the response is given by Eq. I with 
B = sgn mx{by A G^) where 2by = k^ — ky . 

Now let us consider a more generic/isotropic case by 
increasing the tunneling strengths 7 , (3 large enough so 
that the insulating gap closes and a single closed EL in¬ 
side the BZ forms. The gapless semimetallic phases of 
H 3 can be found using the constraints that = 0 or tt 
and (1 =f1 + /3 + 7 —w) = j3 cosky A-'y coskz- For intuition 
we can consider the continuum limit and search for solu¬ 
tions near special points in the BZ. First let us consider 
an expansion around the origin where the constraint is 
I3ky + ^k1 = 2m, i.e., the equation for an ellipse. As¬ 
suming that / 3,7 > 0 to be explicit, this constraint only 
has a solution when m > 0. Now to be concrete take 
/3 = 7 = 2m = 2 so that there is a only a single EL circle 
located in the kx = 0 plane and none at kx = tt. To illus¬ 
trate the nature of the EL let us expand the Hamiltonian 
near it. It is convenient to switch to cylindrical coordi¬ 
nates: {kx, ky, kz) -A {kx, K, 0) where 0 winds around the 
EL and k represents the (signed) radial distance away 
from the EL in the fcyfc^-plane. Using this definition, 
{kx = 0, K = 0,9) lies on the EL and we find the Hamil¬ 
tonian near the EL is 

HFL{k) « Skx(J^+{3/2{6kx)^+2SK)a^ « 6 kxa'^+m{SK)a^ 

( 3 ) 

where the mass function m{SK) = 26 k and nothing de¬ 
pends on 9. Thus, near the Fermi surface we find a fam¬ 
ily of ID Dirac Hamiltonians along the a;-direction with 
masses depending on the radius in k-space away from 
the Fermi-surface {Sk) in the kykz-plane which can be 
positive or negative. This expansion shows that at each 
{ky, kz) we have the Hamiltonian of a massive ID Dirac 
model, and the sign of the mass (and thus topological 
phase) changes as a function of {ky,kz) as one passes 
through the EL. The ID Bloch Hamiltonians Hk ,kSkx) 
along kx (and parameterized by {ky, kz)) will be topolog¬ 
ically trivial on one side of the FL and non-trivial on the 
other side. Which side is topological, and which is triv¬ 
ial depends not only on the nodal lines, but also on the 
secondary weak invariant of the filled bulk bands. This 
bulk weak invariant influences whether the surface states 
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FIG. 1. (a) The polarization for the model in Eq.2 is plotted 
vs the parameter m in the model with /? = 7 = 2. The po¬ 
larization scales approximately linear with m for small values 
of m, but deviates as m is increased, (b) The magnetization 
for the model in Eq.2 with an extra term tpp sin kyl is plotted 
vs various values of m for /3 = 7 = 2. It is linear for small 
values of m, but ends up saturating when m is increased, (c) 
The location of the line node is plotted in the E — ky — kz 
space with kx = 0 for tp = 0.0 and various values of m with 
/? = 7 = 2. The polarization is proportional to the area en¬ 
closed by the FL. (d) The location of the line node is plotted 
in the E — ky — kz space with kx = Q for m = 1 and various 
values of tp with /3 = 7 = 2. The magnetization is propor¬ 
tional to integral of the energy around the FL in momentum 
space. 


in the projected surface BZ are on the interior or exterior 
of the FL. 

As we have seen in the special case above {7 = 0, /3 yf 
0), one quasi-topological EM response determined by the 
geometry of the FLs is the charge polarization. Let us 
consider this more generally. The theory of electric po¬ 
larization for insulators with gapped surfaces [62] was ex¬ 
tended to Chern insulators with gapless edges in Ref. 65, 
and 2D Dirac TSMs in Ref. 42. In the latter two arti¬ 
cles it was shown that insulators and TSMs with gap¬ 
less boundary states can have a well-defined polariza¬ 
tion, though the the connection between the bulk pe¬ 
riodic calculation of the polarization and the boundary 
charge requires care in handling the filling of the bound¬ 
ary modes. For a single FL, the boundary modes can be 
dealt with straightforwardly, but when there are multiple 
FLs the possibility of Z 2 cancellations of overlapping sur¬ 
face states complicates the connection between the bulk 
value of the polarization (derived from Berry phase ar¬ 
guments below) and the surface charge[42]. This can 
be dealt with systematically (as was done for 2D Dirac 
semimetals in Ref. 42), and we resolve the case of two 
FLs in the Supporting Online Material. 

To generally determine the polarization of a LTSM in 
some fixed direction h we consider the family of ID Bloch 


Hamiltonians (fey) parameterized by k_\_, k \\, the com¬ 
ponents of the momentum perpendicular and parallel to 
n. Generically the family iLj:^(A:||) is a set of ID gapped 
Bloch Hamiltonians except when the point {k\\,k±) lies 
on one of the Dirac FLs (which only occupy a set of mea¬ 
sure zero in the 3D BZ). Note that while the FLs in our 
model are planar, our results below apply to non-planar 
cases as well. To calculate the charge polarization we 
first need to calculate the quantity[62] 


e„(fc.) = ^ 


rffciiTr 


A||(fc||,fcj-) 


(4) 


where A\\ is the component of the Berry connection along 
h. Let us first consider the special case where we eval¬ 
uate 0(fcx) at k± = Aa, where is any inversion- 
invariant momentum in the k± plane, i.e., A^ = —A^ 
mod G. Then 0 (A[i) = 0 mod e or | mod e if (fcy) 
is gapped, since this ID Bloch Hamiltonian has inver¬ 
sion symmetry. We then consider a deviation away from 
k± = Aa which is still in the plane normal to n, and such 
that the Hamiltonian -i- 5 fc^(^||) i® gapped. However, 
this ID Bloch Hamiltonian does not have to be inversion 
invariant, and thus it is not immediately obvious how to 
evaluate 0||. However, we can use the following general 
argument to simplify the calculation. Let us evaluate the 
difference in the ID polarizations 


A0|| = 0||(Ao + Slt±) — 0||(Aa) 



J dfcy |Tr |^A||(fc||, Aa + (5fcx) 

'A\\mi,Aa)]} = ^Jjr[F] 


(5) 


where the last expression is a surface integral of the Berry 
curvature 2-form F over the region S bounded by the two 
closed circles located at Aa and Aa -I- Sk± and spanned 
by fey through the cycle of the BZ in the fi direction. 
Since our system has TT symmetry, the only sources of 
Berry curvature are the tt-Aux lines carried by the Dirac 
FLs. Thus, generically A0y = 0 or | depending on the 
parity of the number of Dirac line-nodes enclosed in the 
surface S. In fact, this argument is completely general 
and does not rely on starting at an inversion-invariant 
momentum: A0y = 0 or | only relies on the existence 
of TT symmetry. The ability to start at an inversion- 
invariant momentum informs us that the global constant 
needed to determine the full 0y(fcx) from the knowledge 
of only the A0y(A:x) is either 0 or data which is en¬ 
coded in the secondary weak invariant if Aa is not the 
F-point. For a single FL we see that Bij, and hence, the 
overall charge polarization is simply proportional to the 
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projected area of the FL in the h boundary BZ, i.e., 

eB,j = = f (k^) = (-1)"- ® S n,, 

J±BZ ^ 

( 6 ) 

where flij is the area of the FL projected onto the bound¬ 
ary BZ, S = x(sgnTOx)) and % = ±1 corresponds to the 
FL helicity, i.e., the clockwise/counterclockwise flow of 
the Berry flux along the FL with respect to the normal n. 
This bulk result holds up to the addition of a quantum 
of polarization[62]. Also, changing the secondary weak 
invariant Vij , can change the polarization by a quantum, 
and/or a sign, since it switches the projected area to its 
complement in the surface BZ. For a single FL the effects 
are already taken into account in Eq. 6. For more than 
one FL, the bulk calculation will result in the projected 
areas of all the FLs modulo regions where an even num¬ 
ber of FLs have overlapping projections. As shown in 
the Supplementary Material, when FLs have overlapping 
projected areas the connection between this bulk result 
and the surface charge requires some knowledge of the 
filling of the boundary states. 

When an infinitesimal inversion-symmetry breaking 
mass mx is added to the system, the surface states on 
one side will be filled while the other will be empty. Each 
filled surface state will contribute | charge to the bound¬ 
ary which exactly matches the bound charge required 
from the polarization calculation. We confirm this result 
numerically in Fig. fa where of H^{k) is plotted vs. 
m with the corresponding location of the FL shown in 
Fig. Ic. We choose /3 ,7 = 2 so that there is a single FL 
in the k^ = 0 plane and centered around the origin of the 
BZ. Px should be proportional to the area enclosed by the 
FL given by cosky + cosk^ = 2 — to/ 2. For small values 
of TO, the FL is approximately a circle of radius y/m and 
Px ~ sgnTOxI^. This approximation works well when to 
is small, but underestimates P^ as to is increased. At 
TO = 4, the FL given by cos ky + cos kz = 0 will enclose 
half the area of the BZ. We see that the polarization 
will have the symmetry Px{rn) = ^ — Pxi8 — to) simply 
because when to > 4, the FL is centered around (tt, tt) 
on the boundary BZ. Hence, we will restrict ourselves to 
0 < TO < 4 in Fig. la. 

One corollary of these general arguments is that, while 
it is not forbidden to have just a single closed FL in sys¬ 
tems with TT symmetry, it is forbidden to have only one 
(or an odd number) FLs which traverses a non-trivial cy¬ 
cle of the BZ and meets itself. We can see this because 
calculating any component of the polarization would in¬ 
dicate that the polarization must be opposite on either 
side of the FL, however this is not compatible with the 
periodicity of the BZ and thus must be forbidden. This 
is a 3D line-node generalization of the Fermion doubling 
theorem for Dirac nodes in 2D with TT symmetry. 

Similar to the 2D Dirac TSMs, which have a non¬ 
vanishing orbital magnetization when there is an energy 


difference between the Dirac nodes, LTSMs can also have 
a magnetization that depends on how the band touching 
lines are embedded in energy/momentum space. To pro¬ 
duce this effect in our model we need to change the energy 
along the nodal submanifold, and we can do this, e.g., by 
adding an extra kinetic energy term e(fc)I to P[^{k). Fol¬ 
lowing Refs. 44 and 66, we have 

PfCibc p j3u 

^ y j^lMdbU-\{H,{k) + (7) 

e{k) + E_{k))\dcU-) 

where E-{k), |m_) are the energy and Bloch functions of 
the lower occupied band, and the derivatives are with re¬ 
spect to momentum. This is evaluated in detail in the 
Supplementary Material. The main property that sim- 
plihes this computation is that the Berry curvature is 
composed of just tt-Aux ^-functions due to the TT sym¬ 
metry. We And 




= 



( 8 ) 


where the integration is over the nodal line. Similar to 
the magnetization in the 2D Dirac semimetal, the result¬ 
ing M does not depend on the weak invariant Vij. In a 
generic model, Eq. 20 will include a sum over integrals 
for all distinct FLs. This is a 3D generalization of the re¬ 
sults of Refs. 42 and 45 that relate the magnetization of 
the 2D Dirac semimetal to the energy diAerences between 
the band-touching points. 


We conArm this result numerically by adding an extra 
term tpp sin kyl to H^l^k) and plotting the magnetization 
vs. TO in Fig. lb. Again we have Axed ,d ,7 = 2 so that 
there is only one FL, which has x = +1- The magne¬ 
tization for this case can be evaluated analytically from 
Eq. 20 since the energy only depends on ky. The limits to 
which kz extends for the FL can be calculated using the 
equation for the nodal line {cosky + coskz = 2 — mf2). 
Hence, on the nodal line, ky is a function of kz- The max¬ 
imum value of cos ky = 1 and this means that the maxi¬ 
mum/minimum kz is given by Ekzc = ±cos“^(l — ml2). 
This is valid only when to < 4, while for to > 4, the FL is 
centered around (tt, tt) instead of the origin. The magne¬ 
tization will have the symmetry AT (to) = Mz{8 — to) 
which is why we restrict ourselves to 0 < to < 4 in 
Fig. lb. The magnetization is given by 


M^{m) 

ipp 


sgnTOx 


AttH 


rko 

l-ko 


sin ky dkz ■ 


(9) 


The magnetization is a function of to and does not have 
a simple closed form expression, but has a linear proAle 
in the regime when to is small. For surfaces with low- 
energy modes we can give a microscopic argument for 
the existence of the magnetization. The surface states 
of H^{k) are initially Aat-bands that do not disperse and 
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FIG. 2. Dotted yellow lines represent initial four-fold degen¬ 
erate Fermi-line (SO). Purple and blue solid lines represent 
spin-split Fermi-surfaces {S+, S—) with (a) a majority and 
minority spin Fermi-line induced by certain T-breaking terms 
(b) spin-split Fermi-lines with equal sizes for each spin rem¬ 
iniscent of a Rashba-type splitting from spin-orbit terms in¬ 
duced by strain/inversion breaking. For both panels the gray 
shaded region represents the magnitude of the polarization 
in the a:-direction from the projected areas of the Fermi lines 
after Z 2 overlap cancellation. 


e{k) will impart a dispersion as a function of {ky,kz). In 
general, this will create a bound surface current in the y— 
z plane which is the consequence of a non-vanishing bulk 
magnetization density. There will be similar currents on 
surfaces without low-energy modes, but there is not as 
simple of an interpretation[42]. 

We have now completed our goal of showing that the 
FL EM response is given by Eq. 1. To conclude, we com¬ 
ment on the applicability of our results to real materials. 
The magnetic heterostructure proposed in Refs. 55 and 
56 breaks T explicitly, hence the spins are not degenerate 
and the line nodes occur with just two overlapping bands. 
Thus, this model corresponds precisely to an effectively 
spinless case that has been described throughout this pa¬ 
per and our results can be directly applied. We expect, 
and have confirmed numerically, that this system will 
have a charge polarization. In the case of spin degener¬ 
ate models, which are found, for example, in the carbon 
allotrope materials in Refs. 57 and 58, a further reflec¬ 
tion symmetry is required to stabilize the LTSM arising 
from four overlapping bands as shown in Ref. 67. Eor 
doubly-degenerate bands the charge polarization, being a 
Z 2 quantity, is trivial. However, we can break spin degen¬ 
eracy by including certain T-breaking terms, or inducing 
additional spin orbit terms via strain, with the require¬ 
ment that the FLs are not completely destabilized to a 
gapped, or point-node, phase. If we take two copies of 
our model, one for each spin, then two illustrations of ini¬ 
tially spin-degenerate FLs (in the kx = 0 plane) split by 
two types of spin-dependent terms are shown in Fig. 2. 
In these cases, the polarization Px can be nontrivial and 
is not completely Z 2 canceled. In fact, in both cases, the 
shaded areas correspond to the magnitude of the polar¬ 
ization, assuming a vanishing secondary weak invariant. 


The magnetization, on the other hand, is not a Z 2 quan¬ 
tity and can be non-vanishing even for four-fold degen¬ 
erate FLs. Hence, we expect that these systems would 
exhibit charge polarization when the FLs are spin-split 
via strain or other spin-dependent perturbations. 

We would like to thank P.Y. Chang, V.K. Chua, V. 
Dwivedi, and A. Tiwari for discussions. We acknowledge 
support from ONR YIP Award N00014-15-1-2383. 


[1] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

( 2010 ). 

[2] B. A. Bernevig, Topological Insulators and Topological 
Superconductors (Princeton University Press, 2013). 

[3] F. Haldane, Phys. Rev. Lett. 61 , 2015 (1988). 

[4] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 
78 , 195424 (2008). 

[5] A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. Ludwig, 
Phys. Rev. B 78 , 195125 (2008). 

[6] A. Kitaev, arXiv preprint arXiv:0901.2686 (2009). 

[7] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 

[8] J. C. Y. Teo, L. Fu, and C. L. Kane, Phys. Rev. B 78 , 
045426 (2008). 

[9] L. Fu, Phys. Rev. Lett. 106 , 106802 (2011). 

[10] T. L. Hughes, E. Prodan, and B. A. Bernevig, Phys. 
Rev. B 83 , 245132 (2011). 

[11] A. M. Turner, Y. Zhang, R. S. Mong, and A. Vish- 
wanath. Physical Review B 85, 165120 (2012). 

]12] C. Fang, M. J. Gilbert, and B. A. Bernevig, Phys. Rev. 
B 86, 115112 (2012). 

[13] J. C. Teo and T. L. Hughes, Phys. Rev. Lett. Ill, 047006 
(2013). 

[14] R.-J. Slager, A. Mesaros, V. Juricic, and J. Zaanen, Nat. 
Phys. 9 , 98 (2012). 

[15] W. A. Benalcazar, J. G. Teo, and T. L. Hughes, arXiv 
preprint arXiv:1311.0496 (2013). 

[16] T. Morimoto and A. Furusaki, Phys. Rev. B 88, 125129 
(2013). 

[17] G.-K. Ghiu, H. Yao, and S. Ryu, Phys. Rev. B 88, 075142 
(2013). 

]18] C. Fang, M. J. Gilbert, and B. A. Bernevig, Phys. Rev. 
B 87 , 035119 (2013). 

[19] T. L. Hughes, H. Yao, and X.-L. Qi, arXiv: 1303.1539 
(2013). 

[20] Y. Ueno, A. Yamakage, Y. Tanaka, and M. Sato, Phys. 
Rev. Lett. Ill, 087002 (2013). 

[21] F. Zhang, G. Kane, and E. Mele, Phys. Rev. Lett. Ill, 
056403 (2013). 

[22] P. Jadaun, D. Xiao, Q. Niu, and S. K. Banerjee, Phys. 
Rev. B 88, 085110 (2013). 

[23] D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. 
Cava, and M. Z. Hasan, Nature 452, 970 (2008). 

[24] J. Moore, Nat. Phys. 5, 378 (2009). 

[25] Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, 
A. Bansil, D. Grauer, Y. Hor, R. Cava, et al, Nat. Phys. 
5, 398 (2009). 

[26] T. Zhang, P. Cheng, X. Chen, J.-F. Jia, X. Ma, K. He, 
L. Wang, H. Zhang, X. Dai, Z. Fang, et al, Phys. Rev. 
Lett. 103 , 266803 (2009). 

[27] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 





6 


(2005). 

[28] B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 
314 , 1757 (2006). 

[29] M. Konig, H. Buhmann, L. W. Molenkamp, T. Hughes, 
C.-X. Liu, X.-L. Qi, and S.-C. Zhang, Journal of the 
Physical Society of Japan 77 (2008). 

[30] C.-Z. Chang, J. Zhang, X. Feng, J. Shen, Z. Zhang, 
M. Guo, K. Li, Y. Ou, P. Wei, L.-L. Wang, et ai, Science 
340 , 167 (2013). 

[31] S.-Y. Xu, C. Liu, N. Alidoust, M. Neupane, D. Qian, 
1. Belopolski, J. Denlinger, Y. Wang, H. Lin, L. Wray, 
et al, Nat. Comm. 3 , 1192 (2012). 

[32] Y. Tanaka, Z. Ren, T. Sato, K. Nakayama, S. Souma, 
T. Takahashi, K. Segawa, and Y. Ando, Nat. Phys. 8 , 
800 ( 2012 ). 

[33] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81 , 109 
(2009). 

[34] X. Wan, A. M. Turner, A. Vishwanath, and S. Y. 
Savrasov, Phys. Rev. B 83 , 205101 (2011). 

[35] G. B. Halasz and L. Balents, Phys. Rev. B 85 , 035103 

( 2012 ). 

[36] S. M. Young, S. Zaheer, J. G. Y. Teo, G. L. Kane, E. J. 
Mele, and A. M. Rappe, Phys. Rev. Lett. 108, 140405 
( 2012 ). 

[37] Z. Wang, Y. Sun, X.-Q. Ghen, G. Franchini, G. Xu, 
H. Weng, X. Dai, and Z. Fang, Phys. Rev. B 85 , 195320 
( 2012 ). 

[38] Z. Liu, B. Zhou, Z. Wang, H. Weng, D. Prabhakaran, S.- 
K. Mo, Y. Zhang, Z. Shen, Z. Fang, X. Dai, et al, arXiv 
preprint arXiv: 1310.0391 (2013). 

[39] M. Neupane, S. Xu, R. Sankar, N. Alidoust, G. Bian, 
C. Liu, I. Belopolski, T.-R. Ghang, H.-T. Jeng, H. Lin, 
et al., arXiv preprint arXiv:1309.7892 (2013). 

[40] Z. Wang, H. Weng, Q. Wu, X. Dai, and Z. Fang, Phys. 
Rev. B 88 , 125427 (2013). 

[41] B.-J. Yang and N. Nagaosa, arXiv preprint 
arXiv:1404.0754 (2014). 

[42] S. T. Ramamurthy and T. L. Hughes, Phys. Rev. B 92, 
085105 (2015). 

[43] H. B. Nielsen and M. Ninomiya, Physics Letters B 105 , 
219 (1981). 

[44] J. Shi, G. Vignale, D. Xiao, and Q. Niu, Physical review 
letters 99, 197202 (2007). 

[45] D. Xiao, W. Yao, and Q. Niu, Phys. Rev. Lett. 99, 
236809 (2007). 

[46] A. Zyuzin and A. Burkov, Phys. Rev. B 86 , 115133 

( 2012 ). 

[47] M. Vazifeh and M. Franz, Phys. Rev. Lett. Ill, 027201 
(2013). 

[48] Y. Ghen, S. Wu, and A. Burkov, Phys. Rev. B 88 , 125105 
(2013). 

[49] F. Haldane, arXiv preprint arXiv: 1401.0529 (2014). 

[50] Y. Ran, arXiv preprint arXiv: 1006.5454 (2010). 

[51] T. L. Hughes, H. Yao, and X.-L. Qi, arXiv preprint 
arXiv:1303.1539 (2013). 

[52] S. S. Pershoguba and V. M. Yakovenko, Phys. Rev. B 86 , 
075304 (2012). 

[53] P.-Y. Ghang, S. Matsuura, A. P. Schnyder, and S. Ryu, 
arXiv preprint arXiv: 1406.0232 (2014). 

[54] S. Matsuura, P.-Y. Ghang, A. P. Schnyder, and S. Ryu, 
New Journal of Physics 15 , 065001 (2013). 

[55] A. A. Burkov, M. D. Hook, and L. Balents, Phys. Rev. 
B 84, 235126 (2011). 


[56] M. Phillips and V. Aji, Physical Review B 90, 115111 
(2014). 

[57] K. Mullen, B. Uchoa, and D. T. Glatzhofer, arXiv 
preprint arXiv: 1408.5522 (2014). 

[58] H. Weng, Y. Liang, Q. Xu, Y. Rui, Z. Fang, X. Dai, and 
Y. Kawazoe, arXiv preprint arXiv:1411.2175 (2014). 

[59] Y. Kim, B. J. Wieder, C. L. Kane, and A. M. Rappe, 
Phys. Rev. Lett. 115, 036806 (2015). 

]60] R. Yu, H. Weng, Z. Fang, X. Dai, and X. Hu, Phys. Rev. 
Lett. 115, 036807 (2015). 

[61] J. Zak, Phys. Rev. Lett. 62, 2747 (1989). 

[62] R. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 1651 
(1993). 

[63] X.-L. Qi, T. L. Hughes, and S.-C. Zhang, Phys. Rev. B 
78, 195424 (2008). 

[64] T. L. Hughes, E. Prodan, and B. A. Bernevig, Physical 
Review B 83, 245132 (2011). 

[65] S. Coh and D. Vanderbilt, Phys. Rev. Lett. 102, 107603 
(2009). 

[66] D. Ceresoli, T. Thonhauser, D. Vanderbilt, and R. Resta, 
Phys. Rev. B 74, 024408 (2006). 

[67] C. Fang, Y. Chen, H.-Y. Kee, and L. Fu, arXiv preprint 
arXiv: 1506.03449 (2015). 


Multiple FLs and the polarization 

When we have multiple FLs, the problem of calculating 
the polarization precisely is not quite as simple because 
the boundary charge is decided by the overlap and filling 
of the low energy boundary states that are enclosed by 
the multiple FLs. Despite this, even in the most general 
setting, the polarization can be written down as a signed- 
sum of the various projected areas enclosed by the vari¬ 
ous FLs. As described in the main text, we showed that 
we can perform a simple bulk calculation to determine 
a set of values for these signs. However, a precise sur¬ 
face theorem giving the bound charge associated to the 
polarization change at an interface is meaningful only 
when the occupations of the surface states are specified 
(similar to the complications in Refs. 42 and 65 for the 
polarizations in a Chern insulator or 2D Dirac semi-metal 
respectively). If the boundary occupations are precisely 
known, then one can determine the necessary sign for 
each area contribution that will determine the correct 
surface charge. Hence, the projected areas that deter¬ 
mine the surface charge are decided by the geometry of 
the FLs, but the signs multiplying each area can differ 
from the bulk calculation, and depend explicitly on the 
boundary state occupation. 

The results simplify when there is only one or two FLs 
in the system. In the former case, the surface charge 
is determined (up to a sign decided by the inversion- 
symmetry breaking) by whether the surface states ex¬ 
ist inside or outside the FL. For two (or more) FLs an¬ 
other complication appears due to the possibility of the 
projected areas overlapping in the surface BZ. In these 
cases we can have edge states overlapping, and we expect 
generically that a Z 2 cancelation will occur for the over- 
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FIG. 3. Rules for the modification of xsgnmx for the deter¬ 
mination of the boundary charge for the case of two FLs are 
illustrated. The green shaded areas represent regions where 
edge states exist, and the dark green area represents areas 
where there are overlapping edge states. Case (a) needs a 
reassignment of arrows while case (b) does not. 


lapping states. Now, let us show how we can determine 
the bulk value of the polarization precisely for the case 
of two line nodes. A natural guess for a generalization of 
the polarization formula we have derived in Eq. 6 would 
be 


= (10) 

a 

but this unfortunately does not account for the possi¬ 
ble Z 2 cancelations. To account for this we start off 
by drawing the projected FLs in the appropriate sur¬ 
face BZ perpendicular to the polarization direction. We 
must take care to include arrows indicating the direction 
along which Berry flux is flowing along the FL with re¬ 
spect to the surface normal. The flow is clockwise when 
the product xsgnmx = +1 and counterclockwise for the 
product xsgnmi = — 1 where x corresponds to the FL 
helicity with respect to the normal along the fth direc¬ 
tion. If there are some regions where the projected ar¬ 
eas of the FLs overlap, we have to carefully handle the 
Z 2 cancellation. We assume that any place where two 
FL areas overlap there is a cancellation. We can effec¬ 
tively take into account in our formula after performing 
a simple graphical analysis. First, if the weak invariant 
start off by shading the region around 
(7r,7r), else we leave it unshaded. Then every time we 
cross a FL, we change from shaded to unshaded and vice 
versa. This prescription gives us a unique way of shading 
the entire surface BZ with the projected FLs where alter¬ 
nating regions are shaded. The shaded regions naturally 
represent regions of the surface BZ with stable surface 
states. After we are done with shading, we check if the 


regions which are shaded have an arrow consistently go¬ 
ing clockwise/counterclockwise on its boundary. If they 
do, we sum over the areas of the regions shaded with 
the product xsgn wx for that region coming from the di¬ 
rection of the arrow on the boundary. If the direction 
of arrows is inconsistent, we follow the reassignment of 
the arrows as shown in Fig. 3 and sum over the modified 
areas. 

With more FLs, this prescription does not give us a 
unique answer in regions which have more than two sets 
of edge states overlapping. The sign of the polarization 
arising from these regions depends on the details of how 
the surface states are coupled to give the Z 2 cancelation, 
and hence how the states are occupied. The value of 
the polarization that matches the surface charge is ulti¬ 
mately still a signed sum of the projected areas, but these 
signs can only be determined after the occupation of the 
edge state branches is chosen. All of these issues arise 
due to the Z 2 stability of the edge states, as opposed to 
the Z stable chiral case. We will leave the problem of ex¬ 
haustive treatment of generic FL configurations to future 
work. 


Magnetization in a LTSM 


Let us now calculate the magnetization for our model, 
which will eventually lead us to the generic form for all 
LTSMs. The calculation of the (orbital) magnetization in 
crystalline systems was developed in Refs. 44 and 66, and 
the result of our calculation is essentially an extension of 
the results of the 2D Dirac semimetal shown in Refs. 42 
and 45. To proceed, the adiabatic (Berry) curvatures 
J^xy,J^yz,J^zx for the following generic two-band model 
are calculated: 

H{k) = A{k)a^ + miay + B{k)a^ ( 11 ) 


where mx represents an infinitesimal inversion-breaking 
mass term that must be added to properly calculate the 
magnetization. Note that for the purposes of calculating 
the adiabatic curvatures, the additional e(fc)I term that 
we will add to change the energy of the FL can be ignored 
since its inclusion will not affect the Bloch wavefunctions. 
The adiabatic curvature can be represented by defining 
the unit vector d as 


d{k) 


{A,mx,B) 
-I- -I- 


which yields 


= e^^^dadJbdjd, 


( 12 ) 


(13) 









where di = for i = x,y, z. So for the model in Eq. 11 
we have 


J'Ij ITlx 


diA djB — djAdiB 
{A^ + m\ + i32)3/2 


(14) 


For the case of the semimetal, the limit of mx —>■ 0 
must be taken. Using the identity that lime_>o ^ 2 _^q ,2 = 
x sgnmx 5(a), the curvature can be simplified to 


=7rsgnmi(5(V^^+ S2)- ^^ -. (15) 


If we think about the actual terms A(k) and B{k) from 
the model H^, then we quickly see that the d-function 
only has non-zero support exactly on the line-nodes. 
Generically, when A(k) and B{k) both vanish, then the 
system is gapless (when mx 0), and these gapless re¬ 
gions are the only sources of adiabatic curvature for a sys¬ 
tem with 7~ir symmetry. Thus, in the gapless, semimetal- 
lic limit the only adiabatic curvature in the BZ is local¬ 
ized exactly on the FL, which we know must be the case 
for a model with TT symmetry. 

To finish the magnetization calculation, consider the 
model Hsik) = e(k)I + H^(k) which now has broken T 
and broken T, but preserves TT. The expression for the 
magnetization density in terms of Bloch bands is given 
by[66] 


J -0^lm{d,u.m{k)+E.(k))\d,u.) 

(16) 

where E-(k), |u_) are the energy and Bloch functions of 
the lower occupied band, and the derivatives are with 
respect to momentum. From symmetry, and from the 
fact that the extra kinetic term is proportional to the 
identity matrix, the above expression simplifies to 


Pfabc r j3jL 

= sgnmx^ / j^2e{k)E,,. (17) 

4/i Jbz (27r)'^ 


The expression from Eq. 15 for the curvature can now 
be substituted. Notice that we can do a coordinate 
transformation under the integral from (ka,kb,kc) —>■ 


(ka,A,B) and the Jacobian of the transformation J = 
\diAdjB — djAdiB\ is already sitting in the curva¬ 
ture up to a total sign. Using the property that 
J^S{g{x))f{g{x))\g'{x)\dx = S{u)f {u)du, we can 

rewrite Eq. 17 as 




± sgn mx 


e 

m 


dk^dAdB^ ,xJ(VA^ + B'^) 


( 18 ) 

where the domain of integration has now changed to the 
range of values which A, B take over the BZ and the outer 
signs represent the helicity of the FL, i.e. the sign of the 
Jacobian. We can make a coordinate transformation to 
polar coordinates in A, B ^ r, 9 where we note that r, 9 
could in general depend on 


= isgnmi — 


dk‘^ X rdrd9^ ,x,5(r) 

-—-2 e(k)^ 19 

(iTTj^ r 


which can be simplified by integrating the expressions 
over r, 9. The 5 function localizes the integral to the FL 
and the integral over 9 gives us a factor of 2tt. 


M = ±sgnmx 7 ^ [ e{k)dk (20) 

where we have explicitly indicated that the integration in 
Eq. 20 is over the FL which is equivalent to dR. We note 
that the magnetization results from integrating the en¬ 
ergy of each point on the FL along the line node. Again, 
the ± sign in front of the magnetization tells us the sense 
in which the Berry flux circulates along the string, i.e., 
clockwise or counter-clockwise. This is a simple deriva¬ 
tion of the bulk magnetization in the case of a single line 
node. If there are multiple FLs, contributions to the mag¬ 
netization from each FL using Eq. 20 must be added up, 
but the result is not as complicated as the polarization 
with multiple FLs since the magnetization adds up nor¬ 
mally, not as a Z 2 quantity. It is important to note that 
the connection between the bulk magnetization calcula¬ 
tion and the boundary current can depend on the details 
of how the boundary states are filled similar to what was 
shown in Refs. 42 for 2D Dirac semi-metals. 















